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Quantum spin models have been studied extensively in one and higher dimensions. Furthermore, 
these systems have been doped with holes to study t-J models of SU{2) spin-1/2. Their anyonic 
counterparts can be built from non-Abelian anyons, such as Fibonacci anyons described by SU{2)s 
theories, which are quantum deformations of the SU (2) algebra. Inspired by the physics of SU (2) 
spins, several works have explored ladders of Fibonacci anyons and also one-dimensional (ID) t- 
J models. Here we aim to explore the combined effects of extended dimensionality and doping 
by studying ladders composed of coupled chains of interacting itinerant Fibonacci anyons. We 
show analytically that in the limit of strong rung couplings these models can be mapped onto 
effective ID models. These effective models can either be gapped models of hole pairs, or gapless 
models described by t-J (or modified t-J-V) chains of Fibonacci anyons, whose spectrum exhibits 
a fractionalization into charge and anyon degrees of freedom. The charge degrees of freedom are 
described by the hardcore boson spectra while the anyon sector is given by a chain of localized 
interacting anyons. By using exact diagonalizations for two-leg and three-leg ladders, we show that 
indeed the doped ladders show exactly the same behavior as that of t-J chains. In the strong 
ferromagnetic rung limit, we can obtain a new model that hosts two different kinds of Fibonacci 
particles - which we denote as the heavy r’s and light r’s. These two particle types carry the same 
(non-Abelian) topological charge but different (Abelian) electric charges. Once again, we map the 
two-dimensional ladder onto an effective chain carrying these heavy and light r’s. We perform 
a finite size scaling analysis to show the appearance of gapless modes for certain anyon densities 
whereas a topological gapped phase is suggested for another density regime. 

PACS numbers: 


I. INTRODUCTION 

Quantum statistics is an important aspect of quan¬ 
tum mechanics and it lays down the rules for identify¬ 
ing different classes of particles. In three dimensions, 
the exchange of two particles in a quantum system can 
result in a phase change of either 0 or tt for the wave 
function, leading to bosons and fermions. The scenario 
is very different in two-dimensions (2D) and it is well 
known that such quantum systems give rise to anyonic 
quasiparticle excitations.These anyons are observed 
as excitations, localized disturbances of the ground state 
in systems with topological order. Rich behavior emerges 
especially for non-Abelian anyons, for which the exchange 
of two anyons is described by a unitary matrix. 

Models of interacting non-Abelian anyons draw mo¬ 
tivation from quantum spin models, and in particular 
the Heisenberg model, which has been studied for a wide 
range of lattices including one-dimension (ID) chains and 
ladders of SU ( 2 ) spins. Non-Abelian anyons can be de¬ 
scribed as so-called quantum deformations of the usual 
SU{2) spins. Formally, they are described by SU{2)k 
Chern-Simons theories. These SU{2)j^ theories are ob¬ 
tained by ‘deforming’ the SU{2) algebra so as to retain 
only the first k 1 representations. The most widely 
studied models are the k = 2 and k = 3 theories that 
describe Tsing’ and ‘Fibonacci’ anyons respectively. 

Non-Abelian anyons are expected to be present in 
topologically ordered systems including certain fractional 
quantum Hall states,p-wave superconductors,^^ spin 


models,solid state heterostructures^^”^^ and rotat¬ 
ing Bose-Einstein condensates.^^ In particular, Fibonacci 
anyons occur as quasi holes in the Z 3 Read-Rezayi state 
which might be able to describe the z/ = 12/5 fractional 
quantum Hall state.While several experiments have 
shown evidence for emergent Majorana modes,ex¬ 
perimental evidence of Fibonacci anyons is yet to come. 

In the absence of interactions, non-Abelian anyons 
have an exponentially degenerate ground state manifold. 
On this ground state manifold exchange (braiding) of 
non-Abelian anyons acts in a non-commutative way, by 
bringing about a non-trivial transformation in the degen¬ 
erate manifold of the many-quasiparticle Hilbert space. 
This has generated interest since it can be used for topo¬ 
logical fault-tolerant quantum computation. 

Interactions between non-Abelian anyons can be mod¬ 
eled by generalizations of the spin-1/2 Heisenberg model. 
These models have been studied for chains of Fibonacci 
anyons^^ including nearest neighbour couplings^^ and 
also longer range couplings^^ have been studied. Chains 
of higher spin quasiparticles have also been explored, 
critical phases have been identified and the correspond¬ 
ing conformal field theory (CFT) has been obtained from 
numerical simulations.^^ The effect of disorder for chains 
of Fibonacci anyons has been investigated.^^’^^ 

As a step towards understanding the collective behav¬ 
ior of itinerant non-Abelian anyons, these chains can be 
doped with mobile holes, inspired by the electronic t- 
J model.Electrons confined in ID can exhibit the 
phenomenon of “spin-charge separation”, where the spec- 
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trum can be interpreted in terms of two independent 
pieces, one arising from electric charge without spin, and 
the other from a spinon without any charge.^^ Analo¬ 
gously low-energy effective t-J models have been ana¬ 
lyzed for the case of doped chains of Fibonacci (and 
Ising) anyons. These models exhibit a fractionaliza- 
tion of the spectrum into charge and anyonic degrees of 
freedom,an extension of a phenomenon that exists 
in Luttinger liquids. 

As a step towards two dimensions, anyon models have 
been investigated on chains coupled to form so-called 
quantum ladders of non-Abelian anyons, which provide 
anyonic generalizations of the 2D quantum magnets. 

In this paper, we combine the models mentioned above, 
by studying doped quantum ladders of Fibonacci anyons 
consisting of two or three chains. The structure of the 
rest of the paper is as follows: Section II begins with a 
brief introduction of non-Abelian anyons and introduces 
our model. In section III A, we discuss the strong rung 
coupling limit drawing analogy with the spin ladder mod¬ 
els, listing out all the possible quantum states and the 
phase diagram for the model under study. In section 
IIIB, we list all the phases that arise for weakly cou¬ 
pled rungs. The subsequent sections are dedicated to 
the discussion of all these phases. Section IV discusses 
the hard-core boson models and section V describes the 
golden chain phases. In section VI, we analyse the effec¬ 
tive t — J models for two and three-leg ladders, present¬ 
ing the numerical results for the same. We extend the 
phenomenon of spin-charge separation to (doped) lad¬ 
ders of Fibonacci anyons. In section VII, we introduce a 
new model of heavy and light Fibonacci anyons that is 
obtained for certain density regimes when the rung cou¬ 
plings are FM. We conclude by summarising our results 
in section VIII. 


II. THE MODEL 

A. Introduction to non-Abelian anyons 

In this paper we focus on Fibonacci anyons that are 
described by the SU{2)s theories. SU(2)k Chern-Simons 
theories^^’^^ are so-called quantum deformations of the 
SU{2) algebra. Their degrees of freedom are encoded 
by ‘topological charges’ j, which are generalized angular 
momenta. In contrast to SU{2)^ in SU{2)j^ theories the 
total ‘spin’ j is limited to be j = 0 , • , |. 

Akin to the tensor product of spins, non-Abelian 
anyons can be ‘fused’ according to fusion rules given by 


min{ji+j2,k-ji-j2} 

jl Xj2= V 

j3 = \jl-j2\ 

For example, for the fusion of two anyons with ji ^2 = 
these rules would mean | ( 8 ) | =001 (for k > 2). 
Similarly, for the case of Fibonacci anyons (/c = 3), when 



FIG. 1: (a) Illustration of the standard fusion tree with site 
labels Yi (that can be either r or 1) and bond labels Xi. (b) 
The flat version of the fusion tree, (c) A basis change to a 
different fusion tree using an F-move. 


ii ,2 = 1 the fusion rule reads as 101 = 001. Note 
that this is different from what one would obtain under 
a tensor product of SU{2) spin-1 particles. In the limit 
k oc, however, we recover the SU(2) algebra and the 
rule simply describes the tensor product of two ordinary 
SU{2) spins. 

In the rest of this manuscript, we focus only to the 
Fibonacci theory with /c = 3, unless otherwise stated. 
There the allowed values for the topological charges are 
j = 0, 1, |. But if we look closely at the fusion rules, 

we can make the identification 0 | and 1 ^ Thus, 

the Fibonacci theory has two distinct types of partieles 
which we denote as 1 for the trivial particle with j = 0 
and r for the Fibonacci anyon with j = 1 respectively. 
Using these, the fusion rules read 

101 = 1 

r0l = l0r = T (2) 

r 0 T = 1 0 T. 

We represent a system of N anyons by means of a fu¬ 
sion tree as shown in Fig. 1(a), where the anyonic charges 
of the individual anyons are labelled by Y^. The fusion 
outcome of successive fusion of the anyons are encoded 
by the ‘bond’ labels (links) in the fusion tree, labelled by 
Xi in Fig. 1(a). The constraints on the bond labels due to 
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FIG. 2: Convention for the right handed braid to exchange 
two particles. 


Using these bases the F-matrix is represented as, 


F = 


1 

1 

1 

0-1 0-1/2 
0-1/2 -0-1 


( 4 ) 


fusion rules which must be satisfied at each vertex signif¬ 
icantly reduce the size of the internal Hilbert space. For 
N Fibonacci anyons (r) the Hilbert space grows asymp¬ 
totically as 0^ where 0 = (v/5 +1)/2 is the golden ratio. 
From now on we draw a flat version of the fusion tree, as 
shown in Fig. 1(b). 

To perform an operation on nearest neighbor anyons, 
it is advantageous to change to a different basis in which 
the two-particle fusion outcome is explicit. This is done 
via the F-move shown schematically in Fig. 1(c). The 
F-move depends on the two site labels Yi and and 
the bond labels Xi-i and If at least one of these 

four labels is 1 , then there is only one choice of bond 
labels that satisfies the fusion algebra and the F-move 
is trivial. A non-trivial matrix is obtained only when 
all the four labels are r anyons. Specialising to the case 
where U = the labels for the three bonds 

\xi-i^Xi^Xi-^i) allowed by the fusion rules are 

{|l,r, 1 ), |l,r,r), |r,r, 1 ), |r,l,r), (3) 

which transforms to a new basis |x^_i,after the 
F-move: 

t, l,r), |r, r,r)}. 


where as mentioned above we have a non-trivial 2x2 
submatrix only when also Xi-i = = r. 

Another operation that we need to perform on nearest 
neighbor anyons is that of exchanging (or braiding) them. 
In Fig. 2, we show our convention for a right-handed 
braid. The left-hand braiding is the inverse of the process 
shown here. Under a right-handed braid anyons a and 
b pick up a phase depending on the anyon types, a 
and 6, that are undergoing an exchange and their fusion 
outcome c. Note that whenever a or 6 are 1 , the phase 
is trivial. Non-trivial phases are only obtained are for 
a = b = r: 


^rr ^ g+47ri/5^ ^rr ^ g-37ri/5^ (- 5 ^ 

In order to implement a braid on the standard fusion 
tree, we first have to change basis using an F-move to 
make the fusion outcome of the two anyons explicit, and 
then braid. This process is shown schematically in Fig. 3. 
This is represented by a Braid matrix B acting on the 
bond labels |x^_i,The only non-trivial Braid 
matrix is obtained when both the sites are occupied by 
r anyons. In the basis of Eq. (3) we obtain: 


g4i7r/5 


B = FRF = 


0 

0 

0 

0 


0 0 

g-3i7r/5 Q 

Q g-3i7r/5 

0 0 

0 0 


0 

0 

0 


q_^4i7r/5 _i_ 3i7r/5 

^(e4W5_e-3W5) 


0 

0 

0 


1 /^4i7r/5 g—3i7r/5\ 

03/2 V / 

3i7r/5 _i_ ]_^4i7r/5 

“T 


(6) 


Note that when the two site labels are a 1 and a r, the 
F-moves and the exchange phases are all trivial. The 
Braid matrix is effectively the hopping of the anyon to 
the adjacent site. When both the site labels are 1, the 
Braid matrix is simply given by the identity matrix. 


B. Golden chain model 


In this section we review the so-called ‘golden chains’, 
consisting of ID arrays of localized Fibonacci anyons with 
pairwise interactions between nearest neighbors.In this 
model the Hamiltonian for the magnetic interactions be¬ 
tween anyons is defined in analogy to the Heisenberg ex¬ 
change interaction. We assign an energy — J if the fusion 
outcome of two interacting anyons is trivial. For AFM 
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FIG. 3: Schematic representation of a braid on the fusion tree. 


couplings (J > 0), this favours the fusion outcome of two 
neighbouring anyons to be trivial, while for FM couplings 
(J < 0), the fusion of two anyons is preferred to be r. 
This interaction between nearest neighbor anyons is de¬ 
picted schematically in Fig. 4(a) and is implemented by 
projecting on the identity fusion channel 

/^mag = JVag = -J{FP^F-^). (7) 

where F is the operator corresponding to the F-move 
(see Eq. (4)) and is an operator that projects onto 
the 1 state. In the basis of Eq. (3), the matrix represen¬ 
tation for the (dimensionless) magnetic interaction can 
be written explicitly as: 


h 


mag — 





(8) 


Golden chains with AEM couplings are described by 
the k = 3 restricted solid on solid (RSOS) model which 
is a GET with central charge c = 7/10. Eor ferromagnetic 
couplings the corresponding GET is that of the critical 
3-state Potts model with c = 4/5.^^’^^ 



FIG. 4: (a) Nearest neighbor magnetic interaction of ampli¬ 
tude J. (b) The kinetic hopping (of amplitude t) of an anyon 
to its nearest neighbor vacant site. The blue circles represent 
T anyons while the white circles denote vacant sites or holes. 


C. Itinerant Fibonacci anyons in ID 

To model itinerant anyons we introduce holes, i.e. sites 
with a trivial anyon 1 on some of the sites. The holes and 
r anyons are labelled by different U{1) (electric) charges 
and anyonic (non-Abelian) charges. The r anyons (re¬ 
ferred to simply as ‘anyons’ or ‘r particles’ hereafter) can 
move on the chain which results in an additional kinetic 
energy contribution. This kinetic process is schemati¬ 
cally shown in Fig. 4(b). It involves the hopping of a 
particle, along with its electric and anyonic charge to a 
neighbouring site. 

Our specific model of itinerant anyons is a general¬ 
ization of the electronic t-J model.Assuming a large 
on-site charging energy, we eliminate the possibility of 
doubly occupied sites and allow anyons to only hop to 
empty sites. As in the case of electrons, the low-energy 
effective t-J model allows hopping of the anyons to near¬ 
est neighbor vacant sites and an exchange interaction be¬ 
tween nearest neighbor anyons analogous to the Heisen¬ 
berg interactions explained above in Sec. IIB. The ki¬ 
netic term can then be written as Fkm = —where 
/ikin is the (dimensionless) operator corresponding to the 
nearest neighbor hopping process shown in Fig. 4(b). A 



FIG. 5: Two-leg ladder: (a) Interactions along the leg and 
rung directions, (b) The zig-zag fusion path. The couplings 
have been indicated. 
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FIG. 6: Evaluation of Hamiltonian terms on a two-leg ladder along the leg direction that involve the braid operation. The 
blue circles represent r particles white the white circles are for vacant sites. The green ellipses denote the particles that are 
interacting along the leg direction, (a) the magnetic term (b) the kinetic term. 


t-J chain of itinerant anyons was studied for Ising and 
Fibonacci anyons,revealing a separation of excita¬ 
tions into charge and anyonic excitations, similar to spin- 
charge separation in its electronic counterpart. 


D. Ladders of Fibonacci anyons 

1. Undoped ladders 

Ladders of Fibonacci anyons are formed by coupling 
chains of localized anyons.Anyons interact with their 
nearest neighbors along the leg and rung directions via 
Jieg and Jrung- Thcsc interactions are shown schemati¬ 
cally in Fig. 5(a). As fusion path we choose the zig-zag 
path shown in Fig. 5(b), since it minimizes the effective 
range of interactions on the fusion path. We choose pe¬ 
riodic boundary conditions along the leg direction and 
open boundary conditions along the rungs. 

With this choice of fusion path, nearest neighbor inter¬ 
actions on the rungs are also nearest neighbor along the 
fusion path, while those between anyons on the same leg 
are longer range along the fusion path. Nearest neighbor 
rung interactions can be implemented in exactly the same 
way as for nearest neighbor on a chain (see Sec. IIB). In 
order to evaluate the interactions between r particles on 
the same leg, we have to implement a change of basis, this 
time by braiding them in a clockwise manner until they 
are neighbours along the fusion path. This braiding is 
performed by the unitary braid matrix B (see Eq. (6)). 
Once the r particles are nearest neighbours along the 
path they can interact with the same term as discussed 
above in Sec. IIB. After carrying out the interaction, the 
anyons have to be braided back to their original positions. 

More specifically, for a two-leg ladder, adjacent anyons 
along the leg direction are next nearest neighbours along 
the fusion path (see Fig. 5(b)). Thus, one needs to im¬ 
plement one braid operation. The Hamiltonian for the 


magnetic interactions between nearest neighbor rungs r 
and r + 1 on the upper leg is given by 

(-^mag)r = Aeg^2r-1 (^mag)2r^2r-l 5 ( 9 ) 

and on the lower leg as 

(-^mag)r = Aeg^2r+1 (^mag)2r^2r+l 5 (10) 

where /imag has been defined in Eq. (8) and r labels the 
rungs (so that i = 2r labels the diagonal bonds along the 
path). In Fig. 6(a), we summarize the magnetic inter¬ 
actions between nearest neighbor along the leg direction 
for a two-leg ladder. 

For a three-leg ladder, the leg interactions are longer 
ranged interactions, since adjacent anyons on a leg are 
third neighbours along the fusion path. Thus, nearest 
neighbor leg interactions on a three-leg ladder ladder re¬ 
quire two braids before particles are nearest neighbors on 
the fusion path.^^ One gets for interaction between rungs 
r and r + 1 on the upper leg, 

(^mag)r = Aeg^3r-2^3r-l (hmag)3r^3r-l^3r-2 5 (H) 

on the middle leg, 

(^mag)r = Aeg^3r+1^3r-l (hmag)3r^3r-1^3r+l 5 (12) 

and on the lower leg, 

(-^mag)r = Aeg^3r+2^3r+l (hmag)3r^3r+l^3r+2 5 (13) 

where i = 3r labels the diagonal bonds along the path. 

The full magnetic Hamiltonian on the legs is obtained 
by adding all contributions, i^^lg — YIY ^mag? where W 
is the number of legs. As an implementation detail we 
want to mention that the action of the operator i^^ag can 
generate up to 2^^“^ states for each bond interaction. 
The operator mixes spin labels, thus generates multiple 
images for each two-body interaction. This exponential 
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increase in the number of resulting states leads to denser 
Hamiltonian matrices as one increases the width W of 
the ladder and restricts us from exploring larger system 
sizes. 

For ladders of Fibonacci anyons it was shown that sim¬ 
ilar odd-even effects as seen for SU{2) spins continue 
to exist in the limit of strong AFM rungs.AFM cou¬ 
pled ladders of Fibonacci anyons with even number of 
legs are gapped while those with odd number of legs 
are critical and are described by the same CFT as the 
golden chain. On the other hand, Fibonacci ladders with 
FM rung couplings are quite different from their SU (2) 
counterparts.^^ Fibonacci ladders with a width that is a 
multiple of three are gapped since the rungs form singlets 
{j = 0). Other widths are gapless since the isolated rungs 
effectively form r’s, thereby yielding a gapless chain as 
the effective low-energy model. This is different from 
SU{2) spins, for which an even number form a singlet 
ground state, and where thus all even width ladders are 
gapped. 


2. Doped ladders 


electrostatic 

energy 

> k 


3 


Oi 


© 



la 

-O' 

=1= 


/ V'coul(g') ~ Vrep{q - Qc)^ 



relevant 

subspace 

I I electric charge q 

Z on the rung 

2e* 3e* ^ 


FIG. 7: Schematic energy spectra in the presence of a 
parabolic Coulombic charging energy. One can tune the chem¬ 
ical potential and the repulsion between the particles in order 
to gap out the higher energy sectors. This allows restricting 
the calculations to a low energy subspace. 


In this paper we focus on itinerant ladders of Fibonacci 
anyons. As schematically represented in Fig. 5(a) we de¬ 
note the interaction strengths along the leg direction by 
Jieg and tieg for the magnetic and kinetic terms respec¬ 
tively. Along the perpendicular direction, the couplings 
Aung and trung dcuotc respectively the magnetic and ki¬ 
netic terms. 

The magnetic interactions for the ladder have already 
been described in the section IID 1 and we thus only 
need to discuss the kinetic terms. Once again, the two 
sites on a rung are adjacent on the fusion path and the 
hopping is thus implemented as for the ID t-J chain (see 
section IIC). For a r particle and hole lying on adjacent 
sites on the same leg of the ladder, we need to, once 
again, braid the sites to bring them to adjacent positions 
on the fusion tree. For a two-leg ladder, the kinetic term 
between rungs r and r -f 1 on the upper leg is given by 

(-^kin)'r ~ heg-^2r—1 (^kin) 2 r-^ 2 r— 15 
and on the lower leg as 

(^kin)r = ^leg^ 2r+l(^kin)2r^2r+l: (15) 

where hkin has been defined above and r labels the rungs. 
In Fig. 6(b), we summarize the kinetic process between 
nearest neighbor sites along the leg direction for a two-leg 
ladder. 

For a three-leg ladder, a r particle and hole lying on 
adjacent sites on a leg are third neighbours along the 
fusion path. Thus, nearest neighbor leg interaction on a 
three-leg ladder requires two braids before the r particle 
and the hole are nearest neighbors on the fusion path. 
One gets for kinetic terms between rungs r and r -1- 1 on 


the upper leg, 

(-^kin)’" ~ ^leg-^3r—2-^3r—1 (^kin)3r-^3r—l-^3r—25 (15) 

on the middle leg, 

{^kin)r ~ heg-^3r+l-^3r—1 (^kin) 3 r-^ 3 r—l-^ 3 r+l5 (H") 

and on the lower leg, 

(^kin)r = ^leg^3r+2^3r+l (^kin)3r^3r+1^3r+2, (18) 

where i = 3r labels the diagonal bonds along the path. 

Note that, in the above equations, some of the braids 
may be trivial, in contrast to the case of the magnetic 
interactions. The full kinetic Hamiltonian on the legs is 
obtained by adding all contributions, 

Here, the action of the operator can generate up to 
22 W -2 fQj. each bond interaction. 

We also consider models with an additional rung charg¬ 
ing term 

^Cou\{q) ^ l^ep(^ Qc )t (19) 

where q is the number of anyons on a rung and Qc is de¬ 
termined by the implicit chemical potential. Fig. 7 shows 
the energy profile for a given rung composition on the 
ladder. For the three-leg ladder under consideration, this 
term acts pairwise between all the three possible pairs of 
particles that can exist on the three-site rung. Assuming 
a charging energy which is much larger compared to the 
exchange energy, we can consider the limit Kep ^ oo. 
In this limit we can restrict our calculations to just two 
values of the occupation on the rung, n and n + 1, with 
n ranging from 0 to IF — 1. This reduces the Hilbert 
space and thus allows us to perform simulations of larger 
ladders. 
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FIG. 8: All possible fusion outcomes for a rung of a two-leg 
ladder. The blue circles represent r’s while white circles rep¬ 
resent vacant sites. The U(l) and topological charges for the 
different configurations are denoted in the parenthesis. The 
superscript U{L) refers to the r lying on the upper (lower) 
leg. 


E. Simulation algorithm 

Our numerical simulations have been performed by ex¬ 
act diagonalization using the Lanczos algorithm. We ex¬ 
ploit periodic boundary conditions along the leg direction 
to implement translation symmetry. This allows us to 
block-diagonalize the Hamiltonian into L blocks labeled 
by the total momentum K = 27r^ {m being an integer), 
which reduces the size of the Hilbert space, which is the 
major limiting factor in the simulations, by a factor L. 


III. PHASE DIAGRAMS 
A. Isolated rung limit 

Analogous to standard electronic t-J ladders the 
physics of anyonic t-J ladders can be understood starting 
from the strong rung coupling limit.We thus begin 
by identifying the low-lying states of isolated rungs. 

In Fig. 8, we show the five rung configurations on iso¬ 
lated rungs for a two-leg ladder, and the total U{1) and 
anyonic charges that are possible for these rung config¬ 
urations. When there are two holes on a rung, both of 
these charges are trivially zero. When the rung is occu¬ 
pied by two r particles, the net U{1) charge is 2, however 
the topological charge may be either 1 or r, giving rise 
to two different quantum states |2, r) (named “heavy r”) 
and |2,1) (named “heavy hole” - an empty rung being 
a “light hole”). In the case when there is a single r on 
the rung, it can be either on the upper or on the lower 
leg with charges denoted as (l^,r) and (1^,t) respec¬ 
tively. The corresponding quantum states are respec¬ 
tively and The bonding and anti-bonding 

states (named “light r”) are formed by linear su¬ 

perpositions of the configurations with charges (l^,r) 
and (1 ^,t) given by 

,t)±\1^,t)) . ( 20 ) 




FIG. 9: The various fusion outcomes possible for a doped 
three-leg ladder. The blue circles represent r’s while white 
circles represent vacant sites. The first labels in the paren¬ 
thesis signify the U(l) charge, corresponding to the number 
of r’s present on each rung and the second refers to their fu¬ 
sion outcome. The superscripts U, M, L refer to the positions 
of the T on the different legs (upper, middle, lower) of the 
ladder. 

For a three-leg ladder, many more states are possible, 
as shown in Fig. 9. When there is a single r on the rung 
(^rung = 1): tfi^ “fight t” quantum states formed by 
linear superpositions of the three different positions of 
the r particle are given by 

|l=^,r) = t(|l^,r) + ± \/2|l^,r)) , (21) 

= ^(|l^,r)-|iyr)). (22) 

Depending on the sign of the hopping, one of the states 
|l^,r) acquires the lowest energy. Likewise, when there 
are two r anyons on a rung (Wung = 2) we can form 
quantum states as linear superpositions of the states with 
the same U{1) and topological charges. For J^ng > 0, 
one of the “heavy hole” states 

\/2 + 

± a|2^’^,l)) (23) 

is the ground state, with 

\/^mng ^ ^^rung ~ Jrung 

0^=^ - y. -• (24) 

■^^rung 

If for simplicity, we consider Jmng = Gung, then a = 1. 
For Jrung < 0 either of the “heavy r” states 

|2±,r) = ± V2\2^’^,t)) (25) 

has lowest energy depending on the sign of the hopping 
Gung- When the rung is occupied by three r particles. 
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(a) 


(b) 


FIG. 10: Ground state phase diagram for (a) two-leg ladder and (b) three-leg ladder in the isolated rung limit (with or without 
charging energy). The lowest energy sectors for isolated rungs with different anyon numbers (labelled by A^rung) on a ladders 
of itinerant Fibonacci anyons are indicated. The notations are the same as in Figs. 8 and 9. 


the total U{1) charge on the rung is 3, and the possible 
quantum states are |3,r) (named “super-heavy” r) and 
|3,1) (named “super-heavy hole”) depending on the net 
fusion outcome. 

In the limit of independent rungs and by parametris¬ 
ing the rung couplings as t^ng = cos6> and Jrung = sin6>, 
we have mapped out the parameter space on a unit cir¬ 
cle. In Fig. 10 we show the ground state phase diagram 
for an isolated rung on a two-leg or a three-leg ladder 
for different number A^rung of anyons on the rung. 
Note that these phase diagrams do not depend on the 
charging energy which only gives a constant energy shift 
Kep(^runge* — Qc)^ (depending on the A^rung sector). 


B. Phase diagrams of weakly coupled rungs 

Our goal is now to understand the phase diagram of 
anyonic ladders by starting from the isolated rung cou¬ 
pling limits. We do this by turning on small couplings 
tieg and Jieg between the strongly coupled rungs such 
that |tieg|, I Jiegl << iGungU-^rungl, Order to ensure 
that there is no transition to excited state of the isolated 
rungs. Figures 11 and 12 summarize the phase diagrams 
for two and three-leg ladders. 

Depending on the low-energy states on each rung we 
find six different types of phases: 

• Totally gapped phases (T) appear when there are 


exactly two (for J^ng > 0) or three (for J^ng < 0) 
anyons per rung that fuse into the trivial channel. 
These phases will not be discussed further. 

• Effective golden chains (G^) when there are exactly 
n T anyons on every rung that fuse into a total r. 
An optional ± superscript indicates whether the 
particles are in a bonding (+) or antibonding(—) 
state on a rung. These phases will be discussed in 
Sec. V. 

• Paired phases (P) where two anyons on a rung fuse 
in the trivial channel, forming hard-core bosons. 
These phases will be discussed in Sec. IV 

• Effective t-J chains (G^^) consisting of an effective 
hole that arises from n anyons on a rung fusing in 
the trivial channel and an effective r anyon aris¬ 
ing from m anyons fusing in the r channel. The 
± superscript indicates whether the particles on a 
rung are in a bonding or antibonding state. These 
phases will be discussed in Sec. VI 

• Effective models consisting of two flavors of r 
anyons (P^^) that are formed by fusing m and 
n anyons on a rung respectively. Again the ± su¬ 
perscript indicates whether the particles on a rung 
are in a bonding or antibonding state. 

• A phase separated region PSqs originated from an 
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FIG. 11: Phase diagrams of the two-leg ladder in the strong rung coupling limit, a) is without a rung charging term and b) 
with a large rung charging term Kep- Here the radius denotes the density of anyons. Depending on Filing and coupling several 
phases can be distinguished: a totally gaped phase (T), effective golden chain models (G), effective t-J chains (G), paired 
phases (P), and a phase with two different types of r anyons (D). The legend indicates which rung states are relevant in the 
various phases. See the text for details. 






Gt 

: 12^. r) 


Gs 

: |3,t> 
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Ts 
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FIG. 12: Phase diagrams of the three-leg ladder in the strong rung coupling limit, a) is without a rung charging term and b) 
with a large rung charging term Kep- Here the radius denotes the density of anyons. Depending on Filing and coupling several 
phases can be distinguished: a totally gaped phase (T), effective golden chain models (G), effective t-J chains (G), paired 
phases (P), phase separated phases (PS) and a phase with two different types of r anyons (D). The legend indicates which rung 
states are relevant in the various phases. See the text for details. 


effective t-J chain with dominant attraction be¬ 
tween the effective super-heavy r anyons. 

The effect of a large rung charging energy V^ep is to 
suppress pairing and phase separation in two-leg and 
three-leg ladders. The other phases are unchanged when 
adding this term. We will thus use V^ep = oc to reduce 
the Hilbert space dimension when numerically investigat¬ 
ing the latter phases. 


IV. PAIRING AND EFFECTIVE HARD-CORE 
BOSON MODELS 


effective trivial particle and are then described by mobile 
hard core bosons. In the case of the two-leg ladder this 
phase appears when Jmng > 2|Gung|- Identifying |0,1) 
with an empty site and |2,1) with a hard-core boson we 
end up with an effective hard-core boson (HCB) model, 

^HCB = i +1 + h.c.) + V — 'Ui+i), (26) 

i i 

where b] creates a boson at site i and rii = blbi is the 
boson density. The effective hopping matrix element of 
this hard-core boson model is obtained from second-order 
perturbation theory in tieg to be 


We start our detailed discussion with paired phases 
that arise when two r anyons on a rung fuse into an 


t = —- 


pj-2 


-'rung 


2\t 


rung I 


(27) 
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An effective nearest neighbor attraction between different 
types of holes conies also in second-order and is given by 


V = - 


oj-2 

^heg 


A*ung 2|tr 


(28) 


This is similar to fermionic t-J ladders mapping to a 
Luther-Emery liquid of Cooper pairs. 

In the three leg ladder, similar paired phases described 
by the same hard-core boson model are found when p < 
2/3 and the rung couplings satisfy 


TABLE I: Effective couplings for the various golden chain 
phases as a function of the ladder width W and density p. 


w 

P 

Phase 

J / dleg 

Derivation Eq. 

2 

2 

1 

2 

1 

Gt 

G2 

1 

2 

2 

(A4) 

(A12) 

3 

3 

3 

3 

2 

3 

1 

G? 

G± 

Gs 

3 

8 

11 

8</)2 

1 

(A21) 

(A25) 

(A30) 


30|trung| ^ Aung 


3 


(29) 


The effective hopping matrix element of this hard-core 
boson model is 


VI. EFFECTIVE t-J MODELS 
A. Two-leg ladder 

1. Effective t-J models 


t = — 


oj-2 

^beg 

Ed ’ 


(30) 


where Ed = - 2 V 2 |trung| + and a 

proof is provided in Eq. (A74). An additional effective 
attraction between different types of holes is given by 

V = (31) 


as shown in Eq. (A73). 


V. EFFECTIVE GOLDEN CHAINS 


If all rungs are at the same integer filling n, the effective 
model is either in a trivial gapped phase if the n anyons 
fuse into the trivial channel, or an effective golden chain 
model if they form a total r. We label the latter phases as 
Gn or , where the optional ± index indicates whether 
the anyons are in a bonding state (-F) or antibonding 
state (—) on the rung. 

The phase Gf appears in the two-leg ladder at unit 
filling (outside of the paired phase). With two r particles 
per rung and Aung < 0 we obtain the phase G 2 . 

Analogously, the three leg ladder has effective golden 
chain phases for specific densities p = 1/3, 2/3,1 on the 
ladder. In the case when all the rungs have exactly one r 
particle, the phase Gf is obtained. When there are two 
r anyons per rung and Aung < 0, the phase appears. 
Finally, for three r particles per rung and Aung > 0, we 
obtain the G 3 phase. 

All coupling constants J of the effective Golden chains 
are listed in Table I. 


Doping the p = 1/2 Gf golden chain with (light) holes 
one obtains an effective anyonic t-J chain (G^). The 
magnetic coupling J is the same as for the Gf golden 
chain. 

Increasing the U(l) charge density (p > 1/2), by effec¬ 
tively doping the G^ golden chain with heavy holes, one 
obtains a similar G^ t — J chain for AFM rung couplings 
Aung > 0 with a sign change of the hopping term. 

Coupling constants are summarized in Table II. 


2. Effective model for charge degrees of freedom 

From the mapping of a doped ladder to an effective 
ID t-J chain we expect its spectrum to fractionalize into 
charge and anyon (called also “spin”) degrees of freedom. 
To investigate spin-charge separation in the ladder, we 
first examine the pure 1 /( 1 ) charge spectrum when Aeg = 
0 . 

In the J = 0 limit of an anyoninc t-J chain the itiner¬ 
ant anyons behave like HCBs which can be mapped onto 
a system of spinless fermions. Adding an external flux in 
the ring, the HCB spectrum is therefore given by charge 
excitation parabolas. 


^hcb(p, 0 ext) = — 2 t cos 
j{p) 




. (32) 


where {j{p)} is a set of integers (labelled by the branch 


TABLE II: Effective couplings for the various t-J phases of 
the two-leg ladder as a function of the density p. 


Filling 

Phase 

J / Aeg 

t//leg 

Derivation Eq. 

p<\ 

^01 

1 

2 

1 

2 

(A37) 

p>\ 

^21 

1 

2 

1 

2 

(A42) 
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index p) which determine the continuous momenta, given 
by if = ^ + 5 ) + P(/)ext, P being the density of 

particles in the system. In the J = 0 limit we must be 
careful since the fusion tree labels make the anyons distin¬ 
guishable. Thus, in the absence of magnetic interactions 
the energy levels show a high degree of degeneracy that 
arises due to the built in non-Abelian nature of the Fi¬ 
bonacci anyons. Moving an anyon across the boundary 
cyclically translates the fusion tree labels. All N par¬ 
ticles must be translated over the boundary to be able 
to have the original labeling. This brings about a phase 
shift of (pn = ‘^TrnjN, n being an integer. The charge 
spectrum of the anyonic chain can then be described as 
a union of all HCB spectra for all discrete values of 0^, 
with no external flux: 

KC.r^e = EHCBip,<Pn)- (33) 

The states are labelled by their total momentum Kp ^ = 
Kp ^27r§. 

Our numerical results show that, as expected, the 
charge spectrum of the ladder corresponds exactly to that 
of the effective chain. As an example. Fig. 13 shows the 
Aeg = 0 spectrum of a 2 x 8 ladder with p = 3/4, which 
is in perfect agreement with the spectrum of an effec¬ 
tive J = 0 chain with L = 8 and p = ll2 and with the 
HCB spectrum given by the cosine branches according to 
Eqs. (32) and (33). Note that there is a global shift in 
energy between the ladder and the chain spectra given 

by 

-^shift ~ ^rung-Ns Aung-Nd, (^4) 

where Ns (Nd) are the number of rungs carrying a single 
(two) r (s). 



Momentum 


FIG. 13: Charge spectrum at Jrung = Aung = 1000, tieg = 
1, Aeg = 0,Frep = oo. The solid lines denote the HCB spec¬ 
trum (with an external flux) given by Eq. (32), different colors 
corresponding to the different charge branches (labelled by p 
in Eq. (32). The black circles denote the spectrum of a 2 x 8 
ladder with p = 3/4. The blue crosses correspond to the 
effective chain spectrum for L = 8, y5=l/2 (see Eq. (33)). 


rr 



xj_^^^_I I I I I I_^^^_I I I I ' l 'r- 

100 1000 10000 


leg rung 

FIG. 14: Energy splittings of each of the eight E < 0 en¬ 
ergy levels at A = 0 of the 2x8 ladder with p = 3/4 (sym¬ 
bol/line colors here match the colors of the parabolic branches 
of Fig. 13) as a function of tieg/Aung. The values of the leg 

couplings are tleg — 1^ Aeg ~ 0? Frep ~ OO and trung(— Aung) 

takes different values from 100 to 10000. The fits correspond 
to the expected thing behavior. 


Lastly, we would like to mention that the mapping is 
exact only in the limit when the rung couplings tend to 
infinity. For large, yet finite, rung couplings the energy 
levels of the ladder model in each parabola are split into 
an exponential number of energy levels over a finite en¬ 
ergy range E. This is due to second order processes to 
higher energy states that give rise to a broadening of or¬ 
der E ^ tfgg/trung and tf^g/Aung of the energy levels, as 
shown in Fig. 14. 


3. Numerical comparison between microscopic model and 
ID t-J model 

We next turn on a small Aeg and adiabatically follow 
the splitting of the charge parabolas. Fig. 15(a) zooms 
into the low energy spectrum to show how the highly 
degenerate energy levels are split by a small Aeg = 0 . 1 . 
Fig. 15(b) shows results for a larger coupling Aeg = 1- We 
see that magnetic interactions lift the degeneracy of the 
states with an energy spread proportional to EAeg- This 
is consistent with the behavior of the effective t-J chain 
exhibiting “spin-charge separation”: in Refs. 50,51 we 
showed that the full excitation spectrum of an itinerant 
anyon chain is made up of two independent contributions 
originating from the charge degrees of freedom (described 
in section VIA 2 ) and the anyon degrees of freedom which 
are given by a squeezed (undoped) anyon chain of length 
La = pL where p is the anyon density on a L site t-J 
chain of anyons. 

We now perform a quantitative comparison of the spec¬ 
tra of the anyonic ladder and its corresponding effective 
anyonic chain. As the charge spectra match, we focus on 
the energy difference spectrum (EDS) obtained by sub- 
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Momentum 


FIG. 15: Splitting of the degenerate levels on switching on different values of Jieg on a 2 x 8 ladder with p = 3/4 and 
'/rung = trung = 1000, tieg = 1, Kep = oo. Thc parabolas show the continuous HCB spectrum relevant for Jieg = 0 (see Fig. 13) 
and the red crosses represent the ladder spectrum at (a) Jieg = 0.1 (b) Jieg = 1. 


trading the (supposed) charge excitation component to 
each state. By construction, the EDS then carries the 
information about the anyon degrees of freedom. Note 
that this procedure is only possible at low energy and 
for small enough Jieg, i.e. when a well defined parabolic 
charge branch can be assigned unambiguously to the lev¬ 
els we consider. However, even when a large splitting 
of the energy levels is seen as in Fig. 15(b), we have 
been able to identify exactly the charge excitations cor¬ 
responding to the various levels in the low energy spec¬ 
trum up to an excitation energy of order 5tieg and hence 
obtain the corresponding EDS. The numerical results for 
the EDS on a 2 X 8 ladder with anyon density p = 3/4 
for small Jieg = O.ltieg and intermediate Jieg = tieg are 
shown in Eigs. 16(a) and (b) respectively. We find the 
EDS of the ladder and of the effective chain to be in 
perfect agreement. The perfect mapping of the two-leg 
ladder physics to the physics of the chain hence implies 
straightforwardly that the concept of spin-charge frac- 
tionalization is not strictly ID but also applies to the 
two-leg anyonic ladder, in contrast to the electronic lad¬ 
der analog. 

Note that the EDS subtracted spectrum must not be 


Note that, analogously to the case of the two-leg lad¬ 
der, the mapping to the effective model is exact only in 
the limit of infinite rung couplings. In Eig. 18, we show 
the log-log plot for the broadenings E for each E < 0 
energy levels in the iF = 0 sector as a function of the in¬ 
verse rung couplings for a three-leg ladder. The slope of 
— 1 shows again that, for large but finite rung couplings, 
there are second order effective processes involving higher 
energy states of the rungs. 


confused with the actual energy spectrum of the corre¬ 
sponding squeezed golden chain. In our prescription, sub¬ 
tracting the charge excitations from the full spectrum, 
we get the spectrum corresponding to the anyon degrees 
of freedom as a function of the total momentum of the 
ladder/t-J chain, rather than that of the squeezed golden 
chain. Thus the spectra shown in Eigs. 16(a),(b) are qual¬ 
itatively different from the golden chain spectra. 


B. Three-leg ladder 

All results about the t-J phases are summarized in 
Table III and we provide a short description below. 

1. Density p < 1/3 

Eor p < 1/3, the effective model upon doping the effec¬ 
tive golden chain Gf is again a t-J chain (C^) indepen¬ 
dent of the sign of the couplings. Our numerical spectra 
agree well with the effective model, as shown in Eig. 17 
for a 3x6 ladder with p = 5/18. 


2. Density 1/3 < p < 2/3 and Jrung > 0 

In the density regime 1/3 < p < 2/3, we find an effec¬ 
tive t-J model (C^) upon doping the p =1/3 Gf golden 
chain with heavy holes (increasing the U(l) charge den¬ 
sity) or, equivalently, doping the totally gapped p = 2/3 
T 2 phase with effective light r particles (reducing the 
U(l) charge density). Numerical simulations of three-leg 
ladders in this density regime gives an effective hopping 
that agrees very well with the analytical estimate. 
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Filling 

Aung 

Phase 

J / Aeg 

t/ tleg 

Derivation Eq. 

V|J^e^ 

Derivation Eq. 

p< 1 

any 


3 

8 

1 

(A49) 

— 

— 

1 <p< 1 

> 0 


3 

8 

( 2 +a^) 24 , [(3 + a + 2 V 2 a) (cos 5 ) + l] 

(A55) 

— 

— 

p> i 

> 0 

^23 

1 

1 1 1 q:^ 1 q:^ 1 -1 

2 + 0:2 02 “T 20 ' ' 

(A60) 

2 

(A29) 

p>l 

< 0 

Ct2 

11 

802 

1 1 

20 4 

(A65) 

11 

803 

(A26) 

any 

> 30|trung| 

PSo3 

1 

^(30-2 + 20-3 + 

(A84) 

2 

(A29) 


TABLE III: Effective couplings for the t-J phases of the three-leg ladder. Here Ed = (—0+ |) Aung + V^Lung + , 


3- DGTlSZty p ^ jCLTld Trung ^ 0 

In the strong antiferromagnetic Aung > 0 rung cou¬ 
pling limit, increasing the U(l) charge density starting 
from the gapped chain of heavy holes, introduces 
super-heavy r anyons {i.e. |3,t) states). The system 

is described by a C^s (modified) t-J chain. Unlike in the 
previous simple t-J chains, the super-heavy r’s experi¬ 
ence an effective nearest neighbor attractive potential. 


4 . Density p > 2l3 and Aung < 0 

For ferromagnetic Aung < 0 at a density p > 2/3, we 
map to an effective (modified) t-J chain phase of 
heavy r’s and super-heavy holes. The heavy r’s experi¬ 
ence a nearest neighbor repulsive potential. Our numeri¬ 
cal simulations match very well the analytical estimates. 


5. Phase separation at large Aung > 0 

In the absence of a rung charging energy V^ep an ad¬ 
ditional phase PSqs, that exhibits phase separation, ap¬ 
pears for large Aung > 30|trung|- The physics of the 
interacting light holes and super-heavy r particles is de¬ 
scribed by an extended ID t-J-V model that contains 
the usual couplings of the regular t-J model and, in addi¬ 
tion, the nearest neighbor attractive potential V between 
the super-heavy r particles. The dominant attraction V 
leads to phase separation between an empty and a com¬ 
pletely filled ladder. 


VII. EFFECTIVE MODEL OF HEAVY AND 
LIGHT FIBONACCI ANYONS 

A. The model 

We now discuss a new model which appears for strong 
EM rung couplings on a two-leg ladder with p > 1/2 and 
a large V/ep- A similar effective model also describes the 
three-leg ladder with EM rung couplings and 1/3 < p < 


2/3. When the rung couplings are EM, the fusion of two 
r’s results in a r charge. One thus obtains an effective 
model with two different Fibonacci particles, the heavy 
and light r’s distinguished by their 1/(1) charge. For 
two-leg and three leg ladders, the magnetic interactions 
(similar to the Golden chain) and the potentials between 
the different flavors of r particles are listed in Table IV. 

In addition to the magnetic and potential terms one 
also gets a kinetic process exchanging heavy and light 
r’s on nearest neighbor rungs. This process is shown 
schematically in Fig. 19(a),(b) for the microscopic ladder 
model and the effective chain respectively. Note that, in 
the t-J chain, with holes and r’s, the hopping process 
shown in Fig. 1(c) moves the entire particle along with 
its charges and spin labels. Whereas now, the scenario is 
very different, the spin labels mix with each other as the 
heavy r’s hop over to exchange positions with the light 
r’s. The effective ID model (in the basis of Eq. (3)) for 
the hopping of heavy r’s is described by the Hamiltonian 
^HL given below 

'1 

0 

^HL = / 0 , (35) 

(/)- 3/2 0-1 _ 

where t is the rescaled hopping amplitude. We have 


Width 

ri 

T2 

J / Aeg 

Eq. 


Eq. 

2 

2 

2 

heavy 

heavy 

light 

heavy 

light 

light 

2 

0 

1 

2 

(A12) 

(A16) 

(A4) 

2 

03 

1 

0 

(All) 

(A17) 

3 

3 

3 

heavy 

heavy 

light 

heavy 

light 

light 

11 

802 

5 

80 

3 

8 

(A25) 

(A33) 

(A 21 ) 

11 

803 

5 

80 

(A26) 

(A34) 


TABLE IV: Interactions between heavy and light r’s in two- 
leg and three-leg ladders. The labels n and T 2 indicate the 
type (heavy or light) of the two interacting particles. 
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FIG. 16: Comparison of the energy difference spectra, af¬ 
ter subtracting the charge contribution to the energy of each 
state, of a 2 X 8 ladder, p = 3/4 with that of the effective 
t-J chain L — 8, p = 1/2. The couplings on the ladder are 

Trung — trung — 1000, tleg — I5 Frep — OO and (a) j\eg — 0.1 

and (b) Jieg = 1. 


FIG. 17: Three-leg ladder: The black circles denote the spec¬ 
trum of a 3 X 6 ladder with p = 5/18 while the red crosses are 
for the effective chain with L = 6 and p = 5/6. The couplings 
on the ladder are Jrung = Gung = 1000, tieg = l,14ep = oc. 
(a) Charge spectrum at Jieg = 0. The solid lines denote the 
HCB spectrum, (b) energy difference spectrum for Jieg = 0.1. 


found that, for a two-leg ladder, the effective hopping 
amplitude is 

t = cos(37r/5)tieg. (36) 

A proof of this is provided in Eq. (A47). 

For a three-leg ladder, the same model applies in the 
density regime l/3<p<2/3 (see Fig. 12(b) showing the 
two new low energy states on the rungs) and the hopping 
amplitude is 

t = I [1 + 9 cos(37r/5)] t\eg, (37) 

O 

as shown in Eq. (A71). Note that there exists a sym¬ 
metry between the two kinds of r’s i.e. the number of 
heavy and light r’s in the system can be swapped, leaving 
the physics unchanged. We provide a detailed analytical 
derivation of this Hamiltonian in Appendix A. 


These analytical considerations are found to be in very 
good agreement with our numerical results for Jieg = 0. 
Eig. 20 shows the comparison between the low energy 
spectra of the effective single chain with heavy and light 
taus and the ladder model. Note that for Jieg = 0 
there is an exact E ^ -E symmetry in the spectrum (not 
shown here) that emanates from a symmetry of exchang¬ 
ing heavy and light r's. When there is an odd number of 
both particle types, the momenta are shifted by tt under 
this exchange. 

The ID model allows us to numerically solve larger 
systems with smaller finite size corrections. We, however, 
restrict ourselves to the case Jieg = 0 when there are no 
magnetic interactions between r particles along the leg 
direction but only a small hopping tieg operates between 
the rungs since already this simple model raises several 
open questions. 
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FIG. 18: Energy splittings of all E < 0 energy levels at iF = 0 
of the 3x6 ladder with p = 5/18 as a function of tieg/Gung- 
The values of the couplings are tieg = l,^ieg = 0, Kep = cxd 
and trung(= e/rung) takcs different values from 100 to 10000. 


O Q 


(a) '' 


"flight ^ 


"Hieavy 


- O-.O - 


h 


eg 


-o 


-o 



Tlight Teavy 

^ ^ -O- 


Theavy Tlight 

-• —o- 


FIG. 19: Hopping of heavy (or equivalently light) r’s. (a) 
Microscopic ladder model: the blue circles denote r particles 
and the white circles are vacant sites on the ladder, (b) The 
effective chain. 


B. Single particle dispersion 


FIG. 20: Low energy spectrum of the two-leg ladder with 
strong FM rung couplings compared to the ID heavy and 
light T model. The black circles are for a 2x8 ladder and 
anyon density p = 3/4, while the red crosses are for the ef¬ 
fective chain that it maps to (L = 8,y5 = 1/2). The val¬ 
ues of the couplings are tieg = Ij^/ieg = 0, Kep = co and 

Jrung = -2trung = ”2000. 



We start with the simplest problem of a single light 
(heavy) r moving in a background of L — 1 r’s of heavy 
(light) r. We choose t = — 1 in order to avoid even- 
odd chain length effects (although for L even the sign 
of t is irrelevant). In the spectra shown in Fig. 21 for 
several chain lengths we observe that the dispersion min¬ 
imum is always at iF = 0 and a local minimum appears 
around an incommensurate momentum. The bandwidth 
is about 0.015|tI, independent of whether we consider a 
single heavy or light r. 


C. Critical phase at generic fillings 

We next consider a finite density p of heavy t' s and a 
corresponding filling of 1 — p light r's. Note that due to 
the symmetry between the heavy and light r’s, densities 
p and 1 — p are equivalent. We expect the same behavior 


FIG. 21: Energy dispersions of a single heavy (or light) r 
amidst L — 1 light (or heavy) r’s for even (circles) and odd 
(crosses) length chains, with couplings t = — 1 and J = 0. 

for all densities except for the half-filled case p = 1/2, 
which we shall consider separately in the next section. 
For simplicity we thus choose p = 1/4 since it allows 
us to perform a finite size analysis using three different 
chain lengths L = 12, 16, 20. The corresponding spectra 
are shown in Fig. 22(a). 

One-dimensional gapless systems are often described 
by a CFT and their lowest energy levels are then given 
by 

O'jrn c 

E{L)=eTL+—{--+hL + hR). (38) 

where c is the central charge and Hr are the scal¬ 
ing dimensions of the ‘primary fields’ of the CFT. The 
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FIG. 22: (a) Spectra for the effective chains for heavy and 
light t's at y5 = 1/4 with t = 1, J = 0. Note the momenta 
have been shifted by tt for L odd to make all spectra similar 
and the ground state energy has been subtracted out. The 
lowest energy levels at momenta 0, 27r/L, tt — 27r/L and tt 
have been tagged as A, B, C and D, respectively. The second 
excitations at momenta tt and tt — 27r/L are labeled by E 
and F, respectively, (b) Finite-size scaling analysis of the A, 
B, C, D, E and F energy excitations. Linear (dashed lines) 
and exponential (full lines) hts are shown around K = tt and 
iF = 0, respectively (see text). The scalings of the B, C and 
F gaps are also reported in (a). 


(thermodynamic) ground state energy per site and 
the velocity v are non-universal constants. The finite size 
ground state energy Eq{L) corresponds to = Hr = 0. 

To test the CFT prediction, we performed a finite-size 
scaling analysis of the first few energy gaps vs l/L. As 
shown in Fig. 22(b), we observe that the gaps around 
K = TT show a linear scaling with l/L, suggesting gapless 
modes. This behavior is, in principle, consistent with 
the CFT scaling of Eq. (38). However, at this point, 
we could not identify the CFT that describes our model, 
being limited in Lanczos exact diagonalizations to system 
sizes of less then twenty sites. Using the density matrix 
renormalization group (DMRG) might help to obtain the 



FIG. 23: (a) Spectra of the chains of heavy and light r’s of 
different sizes at p — 1/2 and t = 1, J = 0. The ground state 
energy has been subtracted out and the momenta are shifted 
by TT for odd number of particles of each type so as to get the 
same zero ground state momentum in all cases. The lowest 
energy levels at momenta 0, 27r/L, tt — 27r/L and tt have been 
tagged as B, E, G, A, respectively. The second excitations at 
momenta 0 and tt are labeled by F and G, respectively. The 
third excitation at momentum tt is labeled by D. (b) Finite- 
size scaling analysis of the energy gaps of (a) vs 1/L (see text). 
The scalings of the E and G gaps are also reported in (a). 


central charge, but is left for future studies. 

The energy spectrum around iF = 0 shows a different 
behavior: as shown in Fig. 22(b), the finite size gaps of 
the first excited states at momentum K = t) and K = 
27r/L could be fitted as A(L) = A(oc) + C/Lexp(—L/^), 
where A(oo) 0.04 is a finite energy gap and ^ > 10 
is a correlation length. This suggests that, at density 
p = 1/4, the energy spectrum shows both a gapless mode 
with linear disperson, described by a CFT, and additional 
gapped modes. 


D. Possible topological gapped phase at p = 1/2 

Next we consider the density p = 1/2 where there is an 
equal number of heavy and light r's. We simulated chains 
with lengths L = 14, 16, 18, and 20 and show these spec¬ 
tra in Fig. 23(a), revealing low energy excitations at mo- 
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merit a if = 0 and K = it. Performing a finite size scaling 
analysis on the low lying states using system sizes L rang¬ 
ing from 14 to 20 sites, as shown in Fig. 23(b) we find that 
an exponential form like A(L) = A{oo)-\-C/Lex.p{—L/^) 
provides reasonably good fits of the data. These fits 
suggest that three of the gaps extrapolate to zero and 
the next higher energy excitations extrapolate to a finite 
value A(oo) ^ 0.05. Note however, that the correlation 
lengths extracted from the fits are of the order of the 
system size so that our extrapolations have to be taken 
with caution. However, if correct, our findings would in¬ 
dicate a topological gapped phase with a four-fold degen¬ 
erate ground state, although dimerization is not excluded 
(since ground state momenta are both 0 and tt). In any 
case, we believe half-filling is a special case and very dif¬ 
ferent from the other density regimes we considered. This 
behavior is also notably different from the golden chains 
which are known to be gapless for both FM and AFM 
leg couplings. 


VIII. CONCLUSIONS AND OUTLOOK 

We have studied ladders of itinerant Fibonacci anyons, 
determined their phase diagrams and constructed effec¬ 
tive low energy models. In the limit of strong rung cou¬ 
plings we find several different phases: totally gaped 
phases, paired phases described by hard core bosons, 
golden chain phases, t — J phases that carry r anyons and 
trivial particles and lastly the heavy and light r phase 
that carries two flavors of Fibonacci anyons. We have 
derived effective low-energy ID chain models for all of 
these phases. 

The golden chain phases (G) appear for special com¬ 
mensurate fillings, where each rung of the ladder car¬ 
ries a topological charge r in the same quantum state. 
The ladder can hence be mapped onto a golden chain 
with rescaled couplings. Depending on the sign of the 
rung couplings, we also obtain totally gapped phases (T), 
where all the rungs carry a trivial topological charge. 

We obtain additional paired phases (P) described by 
an effective hard-core boson model, similar to fermionic 
t-J ladders mapping to a Luther-Emery liquid of Cooper 
pairs. This phase is obtained in a two-leg ladder when 
ever anyons on a rung always combine to a totally trivial 
charge. 

The most common phase are described by effective t—J 
models (Gnm)- In the case of a three-leg ladder with 
'>^rung > 30|trung| ^ dominant attraction leads to phase 
separation. All other phases show itinerant behavior. 
The mapping onto a t-J chain shows that spin-charge 
fractionalization, which was initially found for ID chains 
of non-Abelian anyons, also occurs in ladders. 

For FM rung couplings, a new phase occurs in which 
the low-energy model contains Fibonacci anyons carrying 
different electric charges but the same topological charge. 
The corresponding effective model is hence very different. 
This phase (Pmn) is seen when Jrung < 0 and p > 1/2 for 


a two-leg ladder while l/3<p<2/3fora three-leg lad¬ 
der. This model allows magnetic interactions between the 
different flavors of r anyons in addition to an exchange of 
the two particle species, where a heavy and light r swap 
positions. We have analyzed this model in the limit when 
there are no magnetic interactions between the different 
flavors of r particles, since it is very rich in itself leading 
to several open questions, including a potential topolog¬ 
ical gapped phase at half-filling. For other densities, we 
find evidence for gapless modes described by a CFT. At 
this stage, it is unclear what is the central charge and 
the CFT that describes our model. These issues deserve 
further studies. 
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Appendix A: Effective models 

In this Appendix we explain the derivation of the effec¬ 
tive models. We perturbatively derive the nearest neigh¬ 
bor couplings of the effective ID models assuming large 
rung couplings and small leg couplings: | Jrung I, |Uung| ^ 

lAegl? I^legl- 

1. Magnetic and Potential Interactions 

Magnetic interactions will only be present in effective 
models if two neighboring sites each have a total anyonic 
T charge. A rung state with a total anyonic r charge can 
have one (light r), two (heavy r) or three (super heavy 
r) T particles on the rung, as shown in Table V. The 
possible magnetic interactions are marked by pink ar¬ 
rows. Below we investigate all possible combinations for 
two neighboring rungs in two-leg and three-leg ladders, 
respectively and compute the effective interactions from 
first-order perturbation in Jieg. 

a. Two-leg ladder 

For the two-leg ladder the rungs have to be in the light 
T state 

|l±,r) = T(|V,r) ± |l^r)) ^ (Al) 

or in the heavy r state |2,t), where we introduced short 
hand notations \U) and \L) for an anyon on the upper 
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TABLE V: The various rung states classified according to 
their U(l) charge (columns) and topological charge (lines), 
j = 0 (j = 1) corresponds to holes (r particles). The possible 
magnetic interactions (in the present models) are indicated 
by pink arrows and the green arrows represent the possible 
kinetic terms arising in first-order perturbation in tieg. The 
leading kinetic and potential interactions can also arise only 
in second or third-order in tieg, marked by blue and orange 
arrows respectively. 


or lower site of a rung, respectively. We review three 
possible cases below. 

Light r — light r : We first evaluate the matrix ele¬ 
ment ('^l^mlgnl'^) of state 

IV’) = ® "t) 

= \{\UU) ±\LU) ±\UL) + \LL)), (A2) 

where i^mlgn (part of the) Hamiltonian that de¬ 

scribes the magnetic interaction between two anyons once 
they are nearest neighbors along a leg of the ladder (see 
main text). 

The only non-vanishing matrix elements for the mag¬ 
netic interaction in the effective model arise when the 
two anyons are on the same leg of the ladder. Thus, we 
have 

= \{{UU\H^Zl^^\UU) + {LL\Ht%J^LL)) 

= \{UU\H^Z%^\UU). (A3) 

The second step follows since the contributions from mag¬ 
netic interactions on the upper and lower legs of the lad¬ 
der are equal in magnitude. It follows immediately that 
the effective magnetic interaction is half the bare inter¬ 
action: 

J = t Jieg. (A4) 

Heavy r — heavy r : Let us now consider the case 
where both rungs are in the heavy r state. The state j^/;) 
is now defined by 

IV’) = |2,r) (g) |2,t) . (A5) 

We calculate the matrix element (V’l-^fmlgnlV’) for this 
state to obtain 

(V’l^fXniV’) = + { ll \ h ^^ i ^^\ ll )) 

= 2(C/C/K%„|17C/), (A6) 


where the second steps follows since both the terms have 
the same contribution. 

To evaluate the contribution of such a term, we need 
to calculate explicitly the matrix elements for all bond 
labels. The bond labels belong to the set S defined in 
Eq. (3), which we rewrite here for convenience 

S = {|lrl), llrr), jrlr), (A7) 

We denote the initial quantum state (including site and 
bond labels) as: 


l^a) = |V’;Va), (AS) 

where |^ck) G S. The final states after the magnetic pro¬ 
cess are given by : 


= (A9) 


where G S. 

Then, the matrix elements for this process in the basis 
S are given by: 




rrleg = _ /. 
magn 


leg 



,(.-2 0-7/2 

(j)~^ + (j)~^ 

(AlO) 

This effective Hamiltonian is proportional to the 
golden chain Hamiltonian (upto an overall shift) and 
the contribution to the effective coupling on each leg is 
Thus, the effective potential between the two 
heavy r’s is given by 


V = 20-Vieg, (All) 

and the effective magnetic interaction for the entire pro¬ 
cess obtained by combining contributions from both the 
legs (see Eq. (A6)), is thus given by 

J = (A12) 

Light r — heavy r : Finally, we consider the case 
where we have a heavy r and light r on neighboring 
rungs. The state IV^) is now defined as 

IV’) = |1'‘',t) O |2,r). (A13) 

We calculate the matrix element ('^l^mlgnl'^) 
state to obtain 

(V’K^nlV’) = l{{UU\HXn\UU) + {LL\H]^1^,\LL)) 

= {UU\HXn\UU), (A14) 

where the second steps follows since both the terms have 
the same contribution. As before we denote the initial 
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and final states as defined in Eq. (A8) and Eq. (A9) re¬ 
spectively. Then, the matrix elements for this process in 
the basis S are given by: 


rrleg = _ 


eg 


A-l 


A-1 


^-2 _0-5/2 

-0-5/2 0-3 J 

(A15) 

This Hamiltonian matrix is equivalent to the golden 
chain Hamiltonian with the effective coupling 


j = -0-1 j. 


leg 


and an effective potential 

F = 0-ljl, 


eg 


(A16) 

(A17) 


b. Three-leg ladder 


Eor three legs, a rung has a total anyonic r charge if 
it is in one of the (light), |2^,t) (heavy) or |3,r) 

(super heavy) state. Below we consider all possible cases 
for neighboring rungs both with a total anyonic r charge. 

Light r — light r : We first assume that the two 
rungs are both in the state 

|l±,r) = i(|V,r) + |l^r)±^/2|ll'^,r)) 

= l|L)±T|M) + i|[/), (A18) 

where the states jL), |M), \ U) denote the position of the 
anyon on the lower, middle or upper leg respectively, and 
± depends on the sign of trung- As above we calculate 
the matrix element for the state 


11=*=,r) 0 


(A19) 


and obtain 


= T^{UU\Ht%J^UU) 

+ {LL\H^^l^JLL) 

= l{UU\H^^l^^\UU), (A20) 

resulting in an effective coupling 

J = — Tleg- (A21) 

Heavy r — heavy r : Let us now assume the effective 

r’s on the rungs are in the heavy r state 

|2*,r) = i(|V’^,r) + |l^’^r)±V2|V’^r)) 

= i(|t/M) + |ML)±y2|t/L)), (A22) 


where we have used the short hand notation \AB) denot¬ 
ing that the r lies on the leg A and B and the net fusion 
channel outcome being a r is implicit. 

As above we calculate the matrix element for the state 

10) = |2±,r)®|2±,r) (A23) 

= ^{\UM,UM) + \UM,ML) + V2\UM,UL) 

+ \ML,UM} + \ML,ML) + V2\ML,UL) 

+ V2\UL,UM) + V2\UL,ML) + 2\UL,UL)) 

and obtain 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

(A24) 


T({UM,UM\HXn\UM,UM) 

{ML,ML\HXn\ML,ML) 

2{UL,UL\HXn\UL,UL) 

{UM,ML\HX,\UM,ML) 

V2{UM,UL\H]^I^JUM,UL) 

V2{ML,UL\H^^l^^\ML,UL) 

V2{UL,UM\Hl^l^^\UL,UM) 

V2{UL,ML\HXn\UL,ML)). 


The first three terms have two magnetic interactions 
each along the leg direction while all the others have only 
one. Moreover all these magnetic interactions are equal 
in magnitude, the contribution of a single term giving the 
effective coupling J = (j)~‘^J\eg (see Eq. (A12)). Taking 
into account the contributions from all the terms, it fol¬ 
lows immediately that the effectively magnetic coupling 
is 


J = y r Vieg (A25) 

and the effective potential is given by 

V = y 0-Vieg . (A26) 

Super-heavy r — super-heavy r : Let us now as¬ 
sume the states on the rungs with total anyonic charge 
r are both given by |3,r) (super-heavy r’s). Then, we 
calculate the matrix element for the state 

10) = |3,r) 0 |3 ,t) (A27) 

and obtain 

+ {LL\H^^l^,\LL) 

+ (MMK«|g„|MM)). (A28) 
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Magnetic interactions on the upper and middle leg con¬ 
tribute to a potential in the effective Hamiltonian given 
by 

V = -2</.-2 Jieg , (A29) 

while the magnetic process on the lower leg results in an 
effective magnetic coupling given by 


J = Ji 


eg • 


(A30) 


Light r — heavy r : Finally, we consider the case 
when the two r’s on the neighboring rungs are in the 
|l^,r) and the |2^,t) states. The state \ 2 p) is now given 
by 

IV-) = |2±,r)®|l±,r) (A31) 

= ^{\UM) + \ML) + V2\UL)) 

® i(|(7) + |L) + V2|M)) 

= \UM,U) + \UM,L) + V2\UM,M) 

+ \ML,U) + \ML,M} + V2\ML,M) 

+ V2\UL,U) + V2\UL,L) + 2\UL,M). 

As before we calculate the matrix element 




= ^{{UM,U\HXn\UM,U) 

+ {ML,L\HXn\ML,L) 

+ 2{ML,M\HXu\ML,M) 

+ 2{UL,U\H'^l^,\UL,U) 

+ 2{UL,L\HXn\UL,L)). 

(A32) 


The contributions of each of these terms to the effective 
coupling is (—0“^)Jieg (see Eq. (A16)), thus the net ef¬ 
fective coupling for this magnetic process is given by 



(A33) 


and the effective potential is 



(A34) 


2. Kinetic terms 

Whenever the total U(l) charges of the two neighbor¬ 
ing rung states differ by ±1, a hopping process occurs in 
first order in tieg. It is the case for a charge-1 (light) r 
and a charge-0 (light) or charge-2 (heavy) hole, a charge- 
1 (light) T and a charge-2 (heavy) r, a charge-2 (heavy) 
hole and a charge-3 (super heavy) r, a charge-2 (heavy) r 
and a charge-3 (super heavy) hole. All these are marked 
schematically by green arrows in Table. V. Below we in¬ 
vestigate all these possibilities in two leg and three leg 
ladders, respectively. 


a. Two-leg ladder 

Let us first evaluate the effective hopping in a two- 
leg ladder arising when there is a r particle in the state 
|l^,r) and an effective hole in either the empty |e) = 
|0,1) (case 1) or the fully occupied |/) = |2,1) (case 2) 
rung state on adjacent rungs. 

Light hole — light r : First, we need to evaluate the 
matrix element (^i|^Hnl^ 2 ) between the two states 

= = (A35) 

and 

| 5 ' 2 ) = |e) O |1=^, t) = \eU)^eL) ^ 

where is the kinetic part of the ladder Hamiltonian 
living on the legs (see main text). Similar to the case of 
magnetic interactions one gets a factor 1/2 and obtains 
for the effective hopping 

t = “beg. (A37) 

Heavy hole — light r : Next, one evaluates the ma¬ 
trix element between the states 

l^i) = |l^, r) ® I/) = (A38) 

and 

l^2) = |/)®|l^,r) = kL±M, (A39) 

The matrix elements are written explicity as 

(A40) 

The derivation is identical but involves two non-trivial 
braids. The contributions from the hopping on the two 
legs are thus given by 

(/'^I'f^kinl'^/) = e"'‘"/®0"4ieg(lT|/lkin|Tl). 

(A41) 

Adding both the terms we get an overall —1 factor, 
compared to the previous case: 

t = — —tieg. (A42) 

Light r — heavy r : Finally, let us consider two 

nearest neighbor rungs, one carrying a heavy r and the 
other with a light r. Here, the bond labels belong to 
the set S defined in Eq. (A7). We denote the ‘initiaL 
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where / is for final and 1 ^/ 3 ) e S. One needs to compute 
all matrix elements The matrix elements 

for the hopping of a heavy and a light r on the lower leg 
of the ladder are found to be: 

where i stands for initial and |^ck) ^ S. The ‘final’ states 
after the hopping process are given by : 

= \lh2i;^i3), (A44) 


quantum state (including site and bond labels) with a 
light (heavy) r on the first (second) rung as: 

= (A43) 


{lh2i-,lTT\H^^^\li2h;lTT) 
{lh2i-,TTl\H^^^l\li2h; ttI) 
{lh2l-,TlT\H^^^\ll2h;TTT) 
{lh2i-,TTT\H^^^\ll2h;TTT) 


ilege“^''*/®(lT|/lkin|Tl), 

-tlege“^’"*/®0“^(lT|/lkin|Tl), 

-ilege“^’'*/®0“^(lT|ft.kin|Tl), 

{lh2l-,TTT\H^^^\ll2h;TlT) = tlege“^’"*/®(/>“^/^(lT|/lkin|Tl), 
tlege“^”/®(/)“^(lT|/lkin|Tl). 


(A45) 


The matrix elements for hopping on the other leg remain 
the same except for the direction of the braid being re¬ 
versed, i.e. all phase factors ^ 

sum up contributions from both the legs, the phase fac¬ 
tors add up and the entire matrix (in the basis (A7)) is 
written as : 


1 


^eff = t 




0 (/)-V2 

0-1/2 0-2 


where the effective hopping amplitude is 
t = [ cos( 37 r/ 5 ))tieg. 


(A46) 


(A47) 


A characteristic feature of this effective hopping Hamil¬ 
tonian is that it mixes the spin labels. This is remark¬ 
ably different from the simple hopping process between 
a hole and a r. In the text, we have quoted the effective 
Hamiltonian matrix for this process as T^hl described in 
Eq. (35), which is related to i^eff by rescaling it such that 
the 2 x 2 block has eigenvalues 0 and 1. More precisely, 
^HL = with a = 0-1 and 5 = 0-^ (where I is 

the identity matrix). 


b. Three-leg ladder 

For a three-leg ladder, we can derive the relevant ma¬ 
trix elements in a similar way. 


Light hole — light r : Let us start with the sim¬ 
plest case and compute the matrix element (Ti|i 7 ^®^|T 2 ) 
where |Ti) = | 5 )(g)|e), IT 2 ) = |e)(g)| 5 ), \s) = is de¬ 

fined in (A18) and |e) is the empty rung. Using obvious 
notations. 


+ 2{eM\H'^^JMe) 

= {eU\H^^jUe}, (A48) 

Since all the F-moves are trivial because of the holes on 
the rungs and there would be no phase factors due to the 
braidings either, we get the effective hopping: 

t = deg. (A49) 


Heavy hole — light r : The calculation is slightly 
more involved when the effective effective hole state \d) 
involves two r anyons on the rung, 

V 2 + 

where X means a vacant site on the rung at position 
X. The matrix element (Ti|i 7 ^®^|T 2 ) now involves the 
initial and final states. 



22 


l^i) 

1 ^ 2 ) 


|s) ® \d) 

^^^{\uu)+-y^\MU)+m 

\d) ® |s) 

^-^(|t7(/) + ^|I7A/) + |t7l) 


+ a\UM) + V2a\MM) + a\LM) + \UL) + V2\ML) + |LL))(A51) 
+ a\MU) + V2a\MM) + a\ML) + \LU) + V2\LM) + |IL))(A52) 


After expanding both sides one gets, 

= 4(2 I ^2) {(LL\H'^^n\UU) + V2a{MM\H'^^JUU) + 2{UM\H'^^JMU) 

+ {UL\H^^^JLU) + a^MU\H^:^jUM) + V2a{LL\H]:^jMM) 

+ V2a{UU\H^:^JMM) + (M L\H^^^JLM) + {LU\H^:^JUL) 

+ 2(LM|F;;^8|ML) + \/2a(MM|F;-8|LL) + (f7!7|F;^8|LL)). (A53) 


The contributions of the individual terms are as follows: 


{LL\H^Z\UU) = {UU\H^:^\LL) = 0-4ieg(lr|/ikin|rl), 

{MM\Hlt\UU) = {UM\h]:^JMU) = {LU\HlfjUL) = e""/Viieg(lr|/iki„|rl), 
{UL\H^^^\LU) = {MM\H^^\LL) = e-4-/®(/.-'tieg(lr|/iki„|rl), 
{UU\H^^\MM) = e4"/50-4ieg(lT|/iki„|Tl), 


{LL\H^^^\MM) 


{MU\H]:^JUM) 

{ML\H^^\LM) 


{LM\HlfjML) = e-4"/Viieg(lr|/ikin|rl). 


(A54) 


Replacing (A54) into (A53) one gets the effective hop¬ 
ping, 


The matrix element for the kinetic Hamiltonian on the 
legs is now expressed as : 


t 


(lr|hkin|rl) 


1 

(2+ 0^2)20 


(3 T o? T 2'\/2q^) 


47r 

cos ^ I + 1 

5 


t\eg 

(A55) 


Heavy hole — super-heavy r : The third case cor¬ 
responds to the effective hole state \d) defined in (A50) 
and the effective r particle state defined by the fully oc¬ 
cupied rung I/) = |3,r). The initial and final states |Ti) 
and 1 ^ 2 ) are now given by : 


1 ^ 1 ) = M)®l/) 

= ijUf) + a\Mf) + \Lf)), (A56) 

1 ^ 2 ) = \f)®\d) 

= ( 1 /^) + «I/M) + I/I)). (A57) 








(A58) 


The individual contributions of these terms are: 
{fU\H^^JUj} = /-2tieg(lr|hki„|Tl), 

(/I|7f;;^8|i:7) = fleg(lT|hkin|Tl). (A59) 

Adding the contributions of the three terms we get: 


t 


(lr|hkin|rl) 


1 

2 + o2 


1 a‘^ o? 


heg- 


(A60) 


Super-heavy hole — heavy r : The fourth case cor¬ 
responds to a r in the state |2^,r) defined in Eq. (A22). 
The effective hole on the ladder is defined by the fully 
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occupied state |/) = |3,1). We calculate the matrix ele¬ 
ment of the states 

|v|>i) = |2±,r)®|/) 

= ^{\UM,f) + \ML,f)±V2\UL,f)){AQl) 

and 

|5'2) = I/) ® |2^,r) 

= ^{\f,UM) + \f,ML)±V2\f,UL)).{AQ2) 

The matrix elements are given by 

+ {f,ML\H'^^JMLJ) 

+ 2{f,UL\H^^^JUL,f)). (A63) 

The contributions of the individual hopping terms are 

{f,UM\H];^JUM,f) = 0-'e®”/®(lr|/ikin|rl), 
{f,UL\H^^jUL,f) = 

{f,ML\H^^^\ML,f) = r'(lr|/iki„|rl). (A64) 

Using Eqs. A64 in Eq. (A63), we get 

(lT|/lkin|'rl) 


Light r — Heavy r : Einally, we consider two rungs 
with a light r defined by the state \s) in (A18) and a 
heavy r defined by the state \h) given as : 

\h) = ^{\L) + V2\M) + \U)). (A66) 

The initial and final states, formed by the tensor prod¬ 
uct of 1 5 ) and \h) are defined as: 

l^'l) = |s)®|/i) 

= ^(\UU) + V2\MU) + \LU) 

+ V2\UM) +2\MM) + V2\LM) 

+ |f/L) + \/2|Ml) + |Ll)) (A67) 

and, 

1 ^- 2 ) = \h)(S>\s) 

= ^(\UU) + V2\UM) + \UL) 

+ V2\MU)+ V2\MM)+ V2\ML) 

+ |If/) + A/2|lM) + |LL)). (A68) 

The matrix elements corresponding to the hopping along 
the leg on the ladder can be expanded using the expres¬ 
sion of the states to give : 


= ^[{LL\H^^^\UU) + 2{MM\H^^\UU) + 2{UM\H^^^\MU) + {UL\H^^^\LU) 

+ 2{MU\H^^JUM) + 2{LL\H^^^\MM) + 2{UU\H]^^\MM) + 2(ML\H^^JLM) 

+ {LU\H'^^JUL) + 2{LM\H^^JML) + 2(MM\H'^^JLL) + (UUIH^^^JLL)). (A69) 

The contributions of the individual terms are as follows: 

{LL\Hlt\UU) = {UU\h]:^JLL) =t,,^{lT\h^,^\Tl), 

{MM\H^^^\UU) = (FM|i7;-s|Mf7) = (I[/|i7i«s|[/I)=iiege-3-/5(lr|/iki„|rl), 

{UL\H^i\LU) = (Ii|7fi[^®|MM) = (Mi|7fi[^®|iM)=tiege3-/5(lr|/iki„|Tl), 

(MMIH^^JLL) = {LM\HifjML)=t,,^e^^^/^{lT\h^in\Tl). (A70) 

Replacing (A70) into (A69) one gets the effective hopping, 

(lT|/lkin|'rl) 

= ^ [1 + 9cos(37r/5)]tieg. (A71) 
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3. Higher order terms 


the intermediate states thus take the value 


In addition to the above cases, we can also have a ki¬ 
netic and potential terms when the difference in the U{1) 
charges on neighboring rungs is larger than 1 . Such pro¬ 
cess occurs e.g. (i) between a charge -0 (light) hole and 
a charge-2 (heavy) hole (marked by blue arrows in Ta¬ 
ble. V) in the P 2 paired phase or (ii) between a charge-0 
(light) hole and a charge-3 (super heavy) r (marked by 
orange arrows in Table. V) in the PSqs phase. 

In case (i), one needs to hop twice to be able to come 
back to a configuration |Tj) that has the same U{1) 
charges as the initial configuration |T/), so the effec¬ 
tive Hamiltonian (leaving in the relevant subspace) is 
obtained in second-order perturbation, in 


1 ^ y ^rung ^^rung 

Ed = -)Jrung — V2trung - - -• (A76) 

Starting from an initial state |3, r) ( 8 ) |e) a simple hopping 
yields the intermediate state | 2 , 1 )( 8 )| 1 , r) with the matrix 
elements given by for hopping on the upper leg 

tjj = (p~^tieg, (A77) 

on the middle leg, 

tu = (A78) 

and on the lower leg. 




^ Er-Ej 


(A72) 

where the sum is on the intermediate states correspond¬ 
ing to two (light) r’s on neighboring sites. The energy 
denominator is given by = (— 2 trung + 'Aung) and ma¬ 
trix elements for all the intermediate states are (—tieg) 
so that one gets for / = J a potential energy 


t\ = tieg . (A79) 

Subsequently, a second hopping yields the intermediate 
state |1,t) (g) 12, 1) with hopping matrix elements on the 
upper leg 


tl = (A80) 


on the middle leg. 




^beg 

Ed 


and for / 7 ^ J a hopping term 


/ = —- 


2i 


leg 


Ed 


(A73) 


(A74) 


In case (ii), one needs to hop three times to be able 
to come back to a configuration that has the same 
1 /( 1 ) charges as the initial configuration, so the effective 
Hamiltonian is now obtained in third-order perturbation 

in 




^ {E, - Ej){Er - Ej) 


(A75) 

where the intermediate states |T^) and carry a 

heavy hole and a light r on neighboring rungs. The 
energy denominators given by the difference in energy 
between the initial (degenerate with the final) state and 


tli = (f> ^iieg, (A81) 

and on the lower leg, 

^2 (A82) 


Finally to come back to a state carrying the same 1/(1) 
charges on the rungs, the hopping amplitudes for the 
three legs are just the complex conjugate of those de¬ 
scribed in Eqs. (A77)-(A79), thus giving for / = J a 
potential energy 


V = 


{24>- 


-2 


1)^2 


leg 


Ed 


(A83) 


and for / 7 ^ J a hopping term 


^ = f(3r 


20-3 + 


(A84) 
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